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5_( ' Abstract. We present and release photometric redshifts for an uniquely large and deep sample of 522286 objects 

I with i'j^g < 25 in the Canada-France Legacy Survey "Deep Survey" fields DI, D2, D3, and D4, which cover a 

total effective area of 3.2 deg^ We use 3241 spectroscopic redshifts with < z < 5 from the VIMOS VLT Deep 
Survey as a calibration and training set to derive these photometric redshifts. Using the "Le Phare" photometric 
redshift code, we devise a robust calibration method based on an iterative zero-point refinement combined with 
a template optimisation procedure and the application of a Bayesian approach. This method removes systematic 
trends in the photometric redshifts and significantly reduces the fraction of catastrophic errors (by a factor of 
2.3), a significant improvement over traditional methods. We use our unique spectroscopic sample to present a 
detailed assessment of the robustness of the photometric redshift sample. For a sample selected at i'j^g < 24, 
we reach a redshift accuracy of cTAz/ii+z) = 0.037 with rj = 3.7% of catastrophic errors (defined strictly as 
those objects with Az/{1 + z) > 0.15). The reliability of our photometric redshifts is lower for fainter objects: 
we find o-Az/(i+z) = 0.029,0.043 and r; = 1.7%, 5.4% for samples selected at i^g = 17.5 - 22.5 and 22.5 - 24 
respectively. We find that the photometric redshifts of starburst galaxies in our sample are less reliable: although 
these galaxies represent only 18% of the spectroscopic sample they are responsible for 54% of the catastrophic 
errors. An analysis as a function of redshift demonstrates that our photometric redshifts work best in the redshift 
range 0.2 < z < 1.5. We find an excellent agreement between the photometric and the VIMOS- VLT deep survey 
(VVDS) spectroscopic redshift distributions at i'^g < 24 for the CFHTLS-Dl field. Finally, we compare the 
redshift distributions of i' selected galaxies on the four CFHTLS deep fields, showing that cosmic variance is 
already present on fields of 0.7 — 0.9 deg^. These photometric redshifts will be made publicly available from 1st 
may 2006 at http://terapix.iap.fr iLnd http://cencosw.oamp.fr/ 
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1. Introduction 

A key factor in the study of galaxy evolution has been 
our ability to acquire large, deep, well-defined redshift 
samples covering substantial volumes of the Universe. 
Since the photometric redshift measurement relies only 
on the measurement of observed colours (Baum 1962), 
this technique can be an efhcient way to assemble large 
and faint samples of galaxies extending to high redshift. 
Moreover, the photometric redshift method is also the only 
way to estimate redshifts beyond the spectroscopic limit 
IjSawicki et al. I L997, 'Arnout s et 10 999, IBenitez E OOO. 
IFontana et"aO 000. rBolzonell ~et al. E 002'l . 

However, this greatly increased redshift-gathering ca- 
pability comes at at price, namely their much lower ac- 
curacy with respect to spectroscopic measurements. The 
most accurate photometric redshifts with medium band 
filters l|Wolf et al. E 004'l still remains around thirty times 
less accurate than redshifts measured with low resolu- 
tion spectroscopy IjLe Fevre et al. E 004b). Despite this, 
for many studies of the galaxy population, such the galaxy 
luminosity function, the velocity accuracy of photometric 
redshifts is sufficient HWolf et al. E 003;i. 

To first order, photometric redshifts are reliable when 
the Balmer or Lyman continuum breaks can be observed 
between two broad band filters. Conventional optical fil- 
ters from B to the / bands can therefore measure red- 
shifts between 0.2 < z < 1. In addition, near infrared 
data are required to provide robust photometric redshifts 
in the "redshift desert" at z > 1.5 since the Balmer 
break is redshif ted t o A > lOOOOA (|Cimatti et al. B 002. 
IGabasch et"aO 004. iMobasher et"arb 004'l. Beyond z > 
3, reliable photometric redshifts can be estimated using 
deep U 01 B band data, based on the Lyman break visi- 
ble at A > 3600A fe.g. lMadauT L995). 

The reliability of photometric redshifts is also re- 
lated to the photometric redshift method. In the stan- 
dard minimisation method, the most likely red- 
shift and galaxy type are determined by a template- 
fitting procedure, which operates by fitting the ob- 
served photometric data with a reference set of spec- 
tral templates (e.g. IPus chell"\982'l. Since no spectro- 
scopic information is required, this standard method 
can provide redshifts beyond the spectroscopic limit 
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l|Bolzonella et al."b 002'). An alternative approach is to 
use a "training method" which can extract information 
from the spectroscopic sample to estimate the photo- 
metric redshifts. For example, neural network methods 
(e.g. IVanzella et al. E 004') or an empirical reconstruction 
of the redshift-colour relation (e.g. [Connolly et al. [ 1995, 
[Wang et al. ( 1998, ICsabai et arb OOOL However, if the 
training set poorly samples the redshift range these meth- 
ods can become unreliable. As a hybrid approach com- 
bining the advantages of both methods, the standard 
method can be optimised using a spectroscopic sample. 
For instance, the initial template set can be optimised 
l|Budavari et al. E OOO. .Bem'tez E 004') or the spectroscopic 
redshift distribution could be introduced as a "prior" in 
a Bayesian fitting procedure IjBenitez E OOOl. Essentially, 
these techniques use the spectroscopic information to im- 
prove photometric redshift quality. 

Until now, however, the major limiting factor in the 
successful exploitation of photometric redshifts has been 
our uncertain knowledge of just how reliable they actu- 
ally are. Are there systematic trends between spectro- 
scopic and photometric redshifts? What fraction of ob- 
jects have 'catastrophic' errors (difference between photo- 
metric and spectroscopic redshifts largely greater than the 
expected uncertainty)? How the photometric redshift reli- 
ability is correlated with the galaxy spectral type and the 
apparent magnitude? Addressing these issues in a thor- 
ough manner requires both a large, highly uniform pho- 
tometric sample free from systematic errors and a large, 
deep, spectroscopic sample of object selected in the sim- 
plest possible manner. Previous studies aX z > 0.3 have 
been concerned either very deep small surveys or larger 
surveys but with correspondingly shallower areas. In all 
cases, the number of available spectroscopic redshifts has 
been small (typically less than ~ 10'^ objects). The combi- 
nation of the high-throughput VIMOS wide-field spectro- 
graph ULe Fevre et al. E003) and t he MEGACAM survey 
camera at CFHT IjBoulade et al. E 003) makes it possible 
to amass a large, deep, highly uniform photometric and 
spectroscopic samples. 

In this paper we present photometric redshifts 
measured using the Canada-France Hawaii Telescope 
Legacy Survey "Deep Fields" catalogues (CFHTLS, 
http://www.cJht.hawaii.edu/Science/CFHTLS) processed 
at the TERAPIX data reduction centre^ complemented 
with shallower VIMOS VLT Deep Survey multi-colour 
data (i^McCracken et al. E 003. ILe Fevre et~al~b 004a'l . We 
use the current release 'T0003' of the CFHTLS. We 
focus on the deep field CFHTLS-Dl (or VVDS-0226- 
04) for which 11567 faint selected spectra I^g < 
24.0 are available from the VVDS spectroscopic survey 
l|Le Fevre et al. E OOSal and are used here as a training 
sample. We then compute photometric redshifts for all 
the CFHTLS "Deep Fields" Dl, D2, D3, and D4 to 
obtain a large and deep dataset of 522286 objects at 
— 25. Photometric and spectroscopic data are de- 

^ terapix.iap.fr 
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scribed in Section 2. Results derived with the standard 
method are presented in Section 3. We describe in 
Section 4 how the standard method can be cahbrated 
using spectroscopic data. The quahty of these cahbrated 
photometric redshifts is described in Section 5, as a func- 
tion of redshift, apparent magnitude and spectral type. In 
Section 6, we investigate how the combination of different 
bands affects the accuracy of our photometric redshifts. 
We finally present in Section 7 the photometric redshifts 
with i'^B < 25 in the 4 CFHTLS deep fields. More detailed 
scientific studies such as the evolution of the angular cor- 
relation function or of the galaxy luminosity function will 
be deferred to forthcoming articles. 

Throughout the paper, we use a flat lambda cos- 
mology {ilm — 0.3, JIa — 0.7) and we define 
h = iJo/100 km s"^ Mpc"^ Magnitudes are given m 
the AB system. Photometric and spectroscopic redshifts 
are denoted by zp and zs; Az represents zp — zs. 

2. Data description 

2.1. CFHTLS multi-colour data 

The MEGACAM deep multi-colour data described in 
this paper have been acquired as part of the CFHT 
Legacy Survey (CFHTLS) which is currently under- 
way at the 3.6m Canada-France-Hawaii telescope. The 
MEGACAM camera consists of 36 CCDs of 2048x4612 
pixel and covers a field-of-view of 1 deg^ with a res- 
olution of 0.186 arcsecond per pixel. The data covers 
the observed wavelength range 3500A < A < 9400A 
in the u* , g' , r', i', z' filters (Figure We anal- 
yse the four deep CFHTLS fields CFHTLS-Dl (centred 
on 02''25"59" - 04°29'40"), CFHTLS-D2 (10''00"28" + 
02°12'30"), CFHTLS-D3 (14''19™27'' + 52°40'56") and 
CFHTLS-D4 (22''15™3P - 17°43'56"), focusing primar- 
ily on the CFHTLS-Dl field for which we have a large 
spectroscopic sample available from the VIMOS-VLT 
deep survey (VVDS). We use the release 'T0003' of the 
CFHTLS. The data processing of the C FHTLS "deep 
fields" is described in IMcCracken et al~T 2006. in prepa- 
ration). Considerable attention has been devoted in the 
TERAPIX pipeline to produce a photometric calibration 
which is as uniform as possible over all fields. Comparing 
the stellar locus in colour-colour planes in the final four 
stacks indicates the variation in absolute photometric 
zero points field-to-field is less than 0.03 magnitudes 
IjMcCracken et al" ~^006. in preparation). 

After removing the masked area, the effective field-of- 
view is about 0.79, 0.80, 0.83 and 0.77 deg^ for CFHTLS- 
Dl, D2, D3 and D4 respectively. In CFHTLS-Dl, the 
catalogue reaches limiting magnitudes of ^ 26.5, 

- 26.4, t'j^q - 25.0, i'^^ - 25.9 and z^^ - 25.0 
(corresponding to the magnitude limit at which we re- 
cover 50% of simulated stellar sources added to the im- 
ages using our default detection parameters). The data 
in other CFHTLS "Deep Fields" are also extremely deep 
with a limiting magnitude 25.7, 26.2, 26.0 in the 



D2, D3, D4 respectively. A summary table listing the 
exposure times in each band is given on the TERAPIX 
web page (http://terapix.iap.fr/). Apparent magnitudes 
are measured using Kron-like elliptical aperture magni- 
tudes (|Kron 1 1980^. The magnitudes are corrected from 
the galactic extinction estimated object by object from 
dust map image s (|Schlegel et al. |L9 98). We multiply all 
the SExtractor l|Bertin <^r^rnoutsl L996'l flux error esti- 
mates by a factor of 1.5 to compensate for the slight noise 
correlation introduced by image re-sampling during the 
stacking of CFHTLS exposures. 

2.2. VVDS multi-colour data 

In addition to CFHTLS data on the CFHTLS-Dl field, 
we use the shallower images from the VVDS survey 
acquired with the wide-field 12K mosaic camera on the 
CFHT l|Le Fevre et al. B 004a'l. McCracken et al. (2003) 
describe in detail the photometry and the astrometry of 
the VVDS-0226-04 field. The VVDS-0226-04 field covers 
the entire CFHTLS-Dl deep field and reaches the limiting 
magnitudes Bab ^ 26.5, Vab ^ 26.2, Rab ~ 25.9 and 
Iab ~ 25.0 (corresponding to 50% completeness). Near 
infrared data in J and Ks bands are also available over 
160 arcmin^ with the magnitude limits of Jab ^ 24.1 
and Kab ~ 23.8 respectively l|Iovino et al. E 005'). 



2.3. VVDS spectroscopic data 

We use the VVDS spectroscopic data acquired with the 
Visible Multi-Object Spectrograph (VIMOS) installed 
at the ESO-VLT. In this paper, we consider the deep 
spectroscopic sample observed in the VVDS-0226-04 
field (CFHTLS-Dl) and selected according to the 
criterion 17.5 < Iab < 24.0 IjLe Fevre et al. E 005a'). 
This sample comprises 11567 spectra. Four classes 
have been established to represent the quality of each 
spectroscopic redshift measurement, corresponding to 
confidence levels of 55%, 81%, 97% and 99% respectively 
IjLe Fevre et al. E 005a'). Since our goal is to assess the 
quality of the photometric redshifts including the fraction 
of catastrophic failures, we restrict ourselves to the classes 
with a confidence level greater or equal to 97% (class 3 
and 4). In the sub- area in common with the CFHTLS-Dl 
field, the final spectroscopic sample used in this paper 
consists in 2867 galaxies and 364 stars with highly 
reliable redshift measurements. The median redshift is 
about 0.76. The la accuracy of the spectroscopic redshift 
measurements is estimated at 0.0009 from repeated 
VVDS observations. 



2.4. Summary 

To summarise, the multi-colour data on the CFHTLS-Dl 
field consists in two joint u*, g' , r', i' , z' and 5, y, i?, 
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Fig. 1. CFHT transmissions curves normalised to miity. 
The solid lines correspond to the u*, g', r', i', z' 
MEGACAM filter curves; the dotted lines correspond to 
the B, V, R, I CFH12K curves. 

/ datasets over 0.79 deg^, adding also J and K appar- 
ent magnitudes over 160 arcmin^. For each object, the 
photometric redshift is computed using all the available 
bands. These photometric redshifts are calibrated using 
2867 spectroscopic redshifts which have a confidence level 
greater or equal to 97%. As an illustration of our com- 
bined photometric and spectroscopic data set, Figure El 
shows the observed colours (using only CFHTLS filters) 
as a function of the spectroscopic redshifts. Multi-colour 
data in u*g'r'i' z' filters is also available on the other four 
fields. 

3. Photometric redshifts with the standard 
method 

We present in this Section the results obtained with a stan- 
dard method, without training the photometric redshift 
estimate on the spectroscopic sample. 

3.1. The photometric redshift code Le^Phare 

We use the code Le Phare^ (S. Arnouts & O. Ilbert) to 
compute photometric redshifts. The standard method 
is described inl Arnouts et al.T l999. 2002). These photo- 
metric redshifts have been found to agree well with com- 
putations from "Hyperz" HBolzonella et al."b OOO). 

The observed colours are matched with the colours pre- 
dicted from a set of spectral energy distribution (SED). 

^ www.lani.oamp.fr/arnouts/LE_PHARE.html 



Each SED is redshifted in steps of Az — 0.04 and con- 
volved with the filter transmission curves (including in- 
strument efficiency). The opacity of the inter-galactic 
medium IjMadau I L995) is taken into account. The merit 
function is defined as 

/=1 V <bs / 

where F^^^JJ^,z) is the flux predicted for a template T 

at redshift z. F^^^ is the observed flux and <t[^^ the as- 
sociated error. The index / refers to the considered filter 
and Nf is the number of filter. The photometric redshift is 
estimated from the minimisation of varying the three 
free parameters z, T and the normalisation factor A. 

3.2. Template set 

Our primary template set are the four Coleman 
Wu and Weedman (CWW) observed spectra: Ell, 
Sbc, Scd, Irr ( |Coleman, Wu fc Weedmair| L980) 
commonly used to estimate the photometric red- 
shifts (jSawicki et al. >L997, IFernandcz -Soto et al"l 999. 
lArnouts et al. 1 999. Brodwin et al. 2006L We add an 
observed starburst SED from [Kinney et al. [ 1996) to 
make our template sets more representative. These 
templates are linearly extrapolated into ultraviolet 
(A < 2000A) and near-infrared wavelengths using the 
GISSEL synthetic models HBruzual et Chariot E 003L For 
spectral types later than Sbc, we introduce a reddening 
E{B -V) = 0,0.05,0.1,0.15,0.2 which follows the inter- 
stellar extinction law measured in the Small Magellanic 
Cloud UPrevot et al. I L984'). Even if these five templates 
are not completely representative of the variety of ob- 
served spectra, it does reduce the possible degeneracies 
between predicted colours and redshift (^Benitez 2 000L 

3.3. Results based on the standard x^ method 

We first apply the standard method on the CFHTLS- 
Dl data without incorporating any spectroscopic infor- 
mation. Figure 13 shows a comparison between the VVDS 
spectroscopic redshifts and the photometric redshifts at 
i' < 22.5. A clear systematic offset is visible at zs < 0.5. 
We would not expect such a trend to appear for such 
a relatively bright sample in a redshift range where the 
Balmer break is between our u* and the r' filters. Small 
uncertainties in the photometric zero-point calibration or 
an imperfect knowledge of the complete instrument trans- 
mission curve are probably responsible for this trend. 

At fainter magnitudes (top left panel of Figure El 
method a)), we sec there are a large number of galaxies 
with Az > 1, mainly in the redshift range 1.5 < zp < 3. 
Most of these catastrophic errors are caused by mis- 
identification of Lyman and B aimer break features. An 
illustration of this degeneracy is presented in Figure 01 
which demonstrates the importance of near-infrared data 
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Fig. 2. Observed colours as a function of the spectroscopic redshifts (black points). The predicted colours derived from 
our optimised set of templates (see section 4.2) are shown with solid lines: Ell (red), Sbc (magenta), Scd (green), Irr 
(cyan) ( [Coleman, Wu &: Weedman | L980) and starburst (dark blue) ( [Kinney et al. [ 1996) from the top to the bottom, 
respectively. 



to break this degeneracy. An alternative solution is to 
include a relevant information in the redshift probabil- 
ity distribution function (PDFz) using the Bayesian ap- 
proach fe.g. lBenfterb OOO. lMoijasher et al."b 004) in order 
to favour one of the two solutions, as is discussed in Section 
Ol 

This basic comparison shows that blindly trusting the 
accuracy of photometric redshifts is perilous. In the fol- 
lowing, we will improve the photometric redshift quality 
using a spectroscopic training set. 

4. An improved method to compute robust 
photometric redshifts 

As we demonstrated in Section 3.3, spectroscopic redshifts 
are required to calibrate the standard method. In this 



Section, we describe the steps we have followed to calibrate 
the photometric redshift estimate. 

4.1. Systematic offsets 

We first select a control sample of 468 very bright galaxies 
{'i'ab — 21.5) which have spectroscopic redshifts. Using a 
minimisation (equation^ at fixed redshift, we deter- 
mine for each galaxy the corresponding best-fitting CWW 
template. We note in each case f/j,^ the observed flux in 
the filter / and F^^^^ the predicted flux derived from the 
best- fit template. For each filter /, we minimise the sum 

Ngal 2 

E ((^/red-^/bs + ^^)K^bs) 
i'<21.5 
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Fig. 3. Comparison between spectroscopic and photomet- 
ric redshifts determined with the standard method 
(without adding the spectroscopic information) for a 
bright selected sample 17.5 < i'j^g < 22.5. 
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Fig. 4. Example of best-fitted templates on multi-colour 
data for a galaxy at zs = 0.311. The solid black 
points correspond to the apparent magnitudes in the 
u* , B, g' ,V,r' , R,i' , I, z' filters from the left to right re- 
spectively. The solid line corresponds to a template red- 
shifted at zp = 2.97 and the dotted line at zp = 0.24. The 
enclosed panel is the associated Probability Distribution 
Function (PDFz). 



filter 



CWW 
i'AB < 20.5 



CWW 
i'AB < 21.5 



CWW 
i *-C 22.5 



PEGASE 
i'AB < 21.5 



u* 


+0.044 


+0.045 


+0.041 


+0.066 


g' 


-0.080 


-0.080 


-0.079 


-0.087 


r 


-hO.Oll 


-0.006 


-0.012 


-0.002 


i' 


-0.005 


-0.002 


-0.004 


-0.001 


t 

z 


-0.037 


-0.025 


-0.014 


-0.045 


B 


-H0.057 


+0.067 


+0.074 


+0.063 


V 


-0.066 


-0.062 


-0.059 


-0.066 


R 


+0.098 


+0.086 


+0.083 


+0.096 


I 


-0.022 


-0.001 


-0.002 


-0.012 



Table 1. Systematic differences between observed and 
predicted apparent magnitudes. These values are given for 
the set of CWW templates and for different cuts in appar- 
ent magnitudes. We add also the values obtained with the 
synthetic library PEGASE. Throughout the paper, we use 
the values quoted for CWW i'^g < 21.5. 



leaving s-^ as a free parameter. For normal uncertain- 
ties in the flux measurement, the average deviation 
should be 0. Instead, we observe some systematic differ- 
ences which are listed in Tabled Such differences have al- 
ready been noted bv lBrodwin et al. 7 2006) in the Canada- 
France Deep Fields Survey. In our data, these differences 
never exceed 0.1 magnitude and have an average ampli- 
tude of 0.042 magnitude. They depend weakly on the 
magnitude cut adopted to select the bright sub-sample 
(TableP) and are also weakly depending on the set of tem- 
plates (see Table^with the values obtained using the syn- 
thetic library PEGASE [Fioc fc Rocca-VolmeraiigeT lQQT). 
Uncertainties in the calibration of the photometric zero- 
points may create discontinuities not reproduced by the 
templates. The size of these systematic differences are 
compatible with the expected uncertainties in the abso- 
lute zero-point calibration (0.05 magnitudes). 

We then proceed to correct the predicted apparent 
magnitudes from these systematic differences, s-^ is the 
estimated correction that we apply to the apparent mag- 
nitudes in a given filter /. If we repeat a second time 
the procedure of template-fitting after having adjusted the 
zero-points, the best-fit templates may change. We check 
that the process is converging: after three iterations the es- 
timated corrections vary less than 2%. The values listed 
in the Tableware measured after 3 iterations. Differently 
from lBrodwin et al. f 2006). the corrections used to correct 
the apparent magnitudes are obtained only on a bright 
sub-sample {i' < 21.5) after 3 iterations. Since the uncer- 
tainties in the zero-point calibration are not better than 
0.01, we add 0.01 in quadrature to the apparent magni- 
tude errors. 
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Fig. 5. Each panel corresponds to one of the four CWW 
templates (Ell, Sbc, Scd, Irr). The points correspond to 
the flux of each galaxy redshifted to the rest-frame using 
the spectroscopic redshifts. The green dashed lines are the 
initial SEDs and the red solid lines are the optimised SEDs 
which are the output of the procedure described in section 
4.2. 

4.2. Template optimisation 

The apparent magnitude measured in the filter Ag/ / pro- 
vides the rest- frame flux at Ag/ //(H- z^) for a galaxy with 
a spectroscopic redshift Zi. Since all the galaxies are at dif- 
ferent redshifts, we can estimate the rest-frame flux over a 
continuous range of rest-frame wavelengths from the spec- 
troscopic sample. In this way can we optimise our set of 
CWW templates. 

We split the galaxy spectroscopic sample according to 
the best-fit template (4 CWW -I- a starburst template with 
a possible additional extinction). Keeping only the objects 
fitted without additional extinction, we use a sub-sample 
of 309 galaxies to perform the template optimisation. The 
black points in Figure [S] show the rest-frame flux recon- 
structed from observed apparent magnitudes. We observe 
a slight deviation between these points and the initial tem- 
plates (dashed lines), particularly for early spectral type 
galaxies. We sort the rest-frame flux according to their 
wavelengths and bin them by group of 50 points. To pro- 
duce the optimised templates, we connect the median flux 
in each bin (solid lines). When no data are available, we 
keep the extrapolation provided by the initial set of tem- 
plates. We don't optimise the starburst template to keep 
the emission lines in this template. 

The colours predicted for these five main optimised 
templates are displayed as a function of redshift in 
Figure El The observed trend in the colour-redshift rela- 



tion are well reproduced by our set of templates. For r' — i' 
and i' — z' colours, we observe oscillations of the predicted 
colour-redshift relation for the starburst template. These 
oscillations are explained by the contribution of emission 
lines like Hq and OIII to the observed flux. Since we use 
only one starburst template to constrain the possible de- 
generacies in colour-redshift space, we are not covering the 
broad range of possible intensities and line ratios. In par- 
ticular, we do not reproduce some blue observed colours 
{r' - i')AB < 0.1 and {i' - z')ab < (Figure EJ. This 
lack of representativeness which we have adopted to avoid 
degeneracies leads to an accumulation of photometric red- 
shifts in certain peaks of the colour-redshift relation, for 
galaxies with strong emission lines. This is responsible for 
the presence of narrow peaks in the redshift distribution 
for the starburst spectral types. 

Finally, these five main optimised templates are lin- 
early interpolated to produce a total of 62 templates to im- 
prove the sampling of the redshift-colour space and there- 
fore the accuracy of the redshift measurement. 



4.3. Bayesian approach 



The Bayesian approach (|Benftez E OOOl allows us to in- 
troduce a relevant a priori information in the PDFz. 
Following the formalism developed bv lBenftez r 2000'). we 
introduce the prior 



p{z,T\i'j,g) =p{T\i'^g)p{z\T,i'^g) 



(2) 



with p{z\T,i'^g) the redshift distribution and p{T\i'j^g) 
the probability to observe a galaxy with the spectral type 
T. p{z\T,i'j^g) is parametrised as: 



, (3) 



p{z\T,i'j^g) cx z"*exp 
and p{T\i'j^g) as: 



Zot + kjntii'AB ~ 20) 



(4) 



The subscript t denotes the type dependency. Using the 
formalism adopted in IBenftez f 2000'). we recompute the 
values of the free parameters using the VVDS redshift dis- 
tribution. We split the sample according to the four opti- 
mised CWW templates. We adjust the parameters at, zot, 
kmt to maximise the likelihood of "observing" the VVDS 
spectroscopic sample. We use the MINUIT package of the 
CERN library IjJames &: Roos I l995) to perform the max- 
imisation (MIGRAD procedure) and to obtain the corre- 
sponding errors (MINOS procedure). The values of these 
parameters for each type are presented in Table |21 The 
parameters ft and kt are also given in Table [21 for types 
one, two, and three and the fraction of type four is auto- 
matically set to complete the sample. 
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spectral type 



at 



zot 



Ell 

Sbc 
Scd 
Irr 



001 +U.iUM 

0.001_Q J^Qg 

1 /1OS+0.081 
1 c;qq+0 038 

1.000_Q ggg 

1 Q4c;+0"21 



0.452 



+0.015 
0.015 
+0.024 



0.166_o.o23 

211+0015 
n 204+0 014 

"•^'-'^-0.014 



n 1 Q7+U.00Y 

Q ,29+0.013 
0.140ing^6 



0.432 
0.080 
0.312 



+0.047 
0.047 
+0.021 
0.021 
+0.033 
0.033 



0.471 

0.306 
0.127; 



0.043 
+0.098 
0.098 
+0.036 



Table 2. Parameters used for the prior P(z,r|i^^) using the formalism from lBenitez r 2000). These parameters are 
derived from the VVDS spectroscopic sample. 



4.4. Summary 

The photometric redshifts are estimated using the code Le 
Phare" (S. Arnouts & O. Ilbert). We calibrate the stan- 
dard method using the VVDS spectroscopic redshifts: 

— We first adjust iteratively the zero-points of the multi- 
colour catalogue using a bright spectroscopic sample. 

— Then we optimise our primary set of templates using 
the observed flux rest-frame shifted at A/(l -f Zi). 

— Finally we apply a prior based on the VVDS red- 
shift distribution following the Bayesian formalism 
presented in IBenite'zT 2000') . 



5. Results: photometric redshift accuracy 

We now assess the quality of the photometric redshifts ob- 
tained with the calibration method described in Section 4, 
by comparing the spectroscopic and photometric redshift 
samples on the CFHTLS-Dl. 



5.1. Method improvement 

Figure El shows the photometric redshifts versus the spec- 
troscopic redshifts for different steps in the calibration 
method. The systematic trends observed with the stan- 
dard -y? method (top left panel, method a)) are removed 
by the template optimisation and the systematic offset 
corrections (top right panel, method h)). After this step, 
the accuracy reaches <Jaz/{i+zs) = 0.037. Adding a prior 
on the redshift distribution decreases the fraction of catas- 
trophic errors without creating any systematic trends 
(bottom left panel, method c J) . The final fraction of catas- 
trophic errors has decreased by a factor 2.3. In the follow- 
ing, we restrict our analysis to the best method c). This 
comparison shows the essential role of the spectroscopic 
information to build a robust photometric redshift sam- 
ple. 

With our final calibration method, we reach an accu- 
racy crAz/(i+zs) = 0.037. At < 24, we recover 96% 
of the galaxies in the redshift range jAzj < 0.15(1 -I- zs). 
(7^2/(1+23) = 0.037 is similar to the accuracy obtained 
by the COMBO-17 survey with a larger set of medium 
band filters (Wol f et al. E 004'). However, considerations on 



the quality of photometric redshifts derived from statisti- 
cal measurements using the whole sample are not really 
meaningful since such statistics depends on the apparent 
magnitude, the spectral type and the redshift range. We 
investigate these dependencies in the next Section. 



5.2. Dependency on apparent magnitude, spectral type 
and redshift 

Figure [71 shows the comparison between photometric and 
spectroscopic redshifts as a function of apparent magni- 
tude. The fraction of catastrophic errors rj increases by 
a factor of 12 going from 17.5 < i'j^g < 21.5 up to 
23.5 < 

i-AB — 24. The redshift rms increases continuously 
from crA2/(i+2s) ~ 0.028 up to aAz/{i+zs) ~ 0.048. The 
apparent magnitude is therefore a key parameter, which 
is to be expected as the template fitting is less constrained 
for the fainter objects. 

Figure jSl shows the comparison between photometric 
and spectroscopic redshifts as a function of the spectral 
type. We define the spectral type according to the best-fit 
template. The fraction of catastrophic errors rj increases 
by a factor of seven from the elliptical to the starburst 
spectral types. The starburst galaxies represent 18% of 
the spectroscopic sample but 54% of the catastrophic er- 
rors. The accuracy in the redshift measurement is sim- 
ilar for the Ell, Sbc, Scd and Irr spectral types with 
<^Az/{i+zs) = 0.032-0.036 but rises to crA2/(i+2s) = 0.047 
for the starburst galaxies. Such a dependency on the spec- 
tral type is expected since the robustness of the photo- 
metric redshifts relies strongly on the strength of Balmer 
break, which is weaker for later types. In addition, the 
photometric redshift estimate of late spectral type galax- 
ies is affected by the intrinsic dispersion in the properties 
of the emission lines and by the large range in intrinsic 
extinction. 

The photometric redshift reliability also depends on 
the considered redshift range. We quantify the dependency 
on the redshift in Figure and Figure EH showing the 
rms scatter crA2/(i+2s) and the fraction of catastrophic er- 
rors 77 as a function of redshift up to z = 1.5. We split 
the sample into a bright 17.5 < i'j^g < 22.5 and a faint 



22.5 < 



'AB 



AB 

< 24. We choose the limit 



■AB 



= 22.5 since 



www.lani.oamp.fr/arnouts/LE_PHARE.html 



it corresponds to the depth of the shallow VVDS and 
zCOSMOS spectroscopic surveys. The fraction of catas- 
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Fig. 6. Photometric redshifts versus spectroscopic redshifts for the sample 17.5 < i'j^^ < 24. Each panel corresponds 
to an additional step in the calibration method with: method a) the standard method ; method b) adding the 
templates optimisation and the corrections of the systematic offsets ; method c) our best method using the Bayesian 
approach, the templates optimisation and the corrections of systematic offsets. The solid line corresponds to zp = zs. 
The dotted lines are for zp = zs zt 0.15(1 + zs). We quote as catastrophic errors the fraction of galaxies with 
\zs — zp\/{l + zs) > 0.15 and the accuracy (Jaz/{i+zs)- The open symbols correspond to galaxies with a second peak 
detected in the PDFz (probability threshold at 5%). 



trophic errors increases dramatically only at z < 0.2. 
At 0.2 < z < 1.5, the accuracy is always better than 
0.045(l + zs)/0.55(l + zs) for the bright and faint samples 
respectively. The fraction of catastrophic errors rj remains 
always less than ~ 4%/14% for the bright and faint sam- 
ple respectively (Figure ITUI) . We observe a degeneracy for 
zs < 0.4 and 1.5 < zp < 3 faint galaxies (bottom left panel 



of Figure . The origin of this degeneracy is a mismatch 
between the Balmer break and the intergalactic Lyman- 
alpha forest depression at A < 1216A. 70% of the galaxies 
at 1.5 < zp < 3 are in fact at zs < 0.4 which prevents the 
use of this spectral range. At zs > 3, the Lyman Break is 
observed between the u* band and the g' bands, allowing 
a reliable photometric redshift estimate for Lyman Break 
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Fig. 9. Accuracy of the photometric redshifts as a func- 
tion of redshift. Catastrophic errors are removed from the 
sample. Only bins with more than ten objects are shown. 

galaxies. We recover 6 of the 8 galaxies at zs > 3 (bottom 
left panel of FigureEJ ■ Even if the quality of the photomet- 
ric redshifts appears good, we point out that we are using 
only spectroscopic redshifts with the highest confidence 
level which is a specific population easier to isolate both 
in photometry and in spectroscopy since they have a sig- 
nificant Lyman break IjLe Fevre et al. E 005b^. Moreover, 
we have tuned the calibration method to be efficient at 
2 < 1.5 using a prior on the redshift distribution (see 
Section 4.3) and without allowing galaxies to be brighter 
than Mbab = -24 Ijllbert et al. B OPS'!. We conclude that 
the most appropriate redshift range for forthcoming sci- 
entific analysis on the complete population of galaxies is 
0.2 <zp< 1.5. 

5.3. Error analysis 

We investigate here the reliability of the error associated 
to the photometric redshift estimate. 

The redshift Probability Distribution Function (see 
lArnouts et al. E 002) is directly derived from the dis- 
tribution 

PDFz = B exp (^-^) , (5) 

with B a normalisation factor. Le_Phare (S. Arnouts & O. 
Ilbert) produces the PDFz for each object. A second red- 
shift solution is likely when a second peak is detected in 
the PDFz above a given threshold. An example of galaxy 
with the good redshift solution enclosed in the second peak 
of the PDFz is shown in Figure 01 The galaxies with a 



Fig. 10. Fraction of catastrophic errors rj per redshift 
bins. The catastrophic errors are defined as galaxies with 
\zs — zp\/{\ -\- zs) > 0.15. The fraction is measured as a 
function of the photometric redshift (top panel) and of 
the spectroscopic redshift (bottom panel). The top panel 
shows the level of contamination, i.e. the fraction of wrong 
redshifts in a given photometric redshift slice. The bottom 
panel shows the level of incompleteness, i.e. the fraction 
of redshifts not recovered in a given photometric redshift 
slice. Only bins with more than ten objects are shown. 



second peak in the PDFz are flagged with open circles in 
Figure El Figured Figure (HI and make up for a large frac- 
tion of the catastrophic errors. We find that the fraction of 
catastrophic errors increases dramatically in those cases: 
when a second peak is detected with a probability greater 
than 5% the fraction of catastrophic errors increases to 
77 = 42%. Removing these galaxies from the sample could 
be useful to select the most robust sub-sample. 

The error bars on the photometric redshifts are given 
by X^{z) = xliin + ^X^- ^X^ = 1 and Ax^ = 9 are 
used to compute the error bars at la and 3cr respectively. 
Fiefure ITTl shows the estimated error bars at 3cr in narrow 
bins of apparent magnitudes. The size of the error bar 
increases towards faint apparent magnitudes, in a consis- 
tent way with the Az rms. We find that 67% and 90% of 
the spectroscopic redshifts are well located in the la and 
3cr error bars respectively. We note that these values re- 
main lower than the theoretical values since photometric 
uncertainties (such as blending or the presence of bright 
neighbours) or the suitability of our template set are not 
taken into account in the computation of the PDFz. We 
conclude that our la error bars are an accurate represen- 
tation of the photometric redshift error, and there can be 
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Fig. 7. Same as Figure with the final cahbration method c). Each panel corresponds to a different selection in 
apparent magnitude. 



a useful way to assess their reliability beyond the spectro- 
scopic limit. 

5.4. Comparison between photometric and 
spectroscopic redshift distributions 

In order to calculate the galaxy redshift distribution, we 
first need to remove the stars from the sample. We use 
the half-light radius ri/2, a morphological criterion mea- 
sured provided by SExtractor IjBertin fc ArnoutsT L996). 
From the spectroscopic sample, we find that 95% of the 
stars have ri/2 < 2.7. Since 16% of the galaxies have also 



ri/2 < 2.7, we combine this morphological criterion with 
a colour criterion. For each object, we compute simulta- 
neously the for the galaxy library and the for the 
star library l|Pickles I L998V If the conditions X^ ~ Xs > 
and ri/2 < 2.7 are satisfied simultaneously, the object is 
flagged as a star. Applying these criteria on the spectro- 
scopic sample, we recover 79% of the stars and only 0.77% 
galaxies are misclassified as stars. The remaining 21% of 
stars are misclassified as galaxies and 69% of these are in 
the redshift range zp < 0.2. 

Since we use the spectroscopic redshift distribution as 
a prior (see Section 4.3) a critical point is to check at which 
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Fig. 8. Same as Figure with our final calibration method c). The sample is selected at 17.5 < < 24. Each panel 
corresponds to a different selection in spectral type defined according to the best-fitting template. 



level the photometric redshift distribution depends on the 
prior. We compare the redshift distributions obtained us- 
ing the prior (weighted solid lines) and without (dashed 
lines) in Figure [T^ The prior has no impact on the global 
shape of the redshift distribution. We see significant dif- 
ferences in the redshift distributions at z > 1.5, where 
the prior efficiently removes the catastrophic failures at 
1.5 < zp <3. 

For the i'j^g < 23 and the i'j^g < 24 selected samples, 
we compare in Figure [T^ the photometric and the VVDS 
spectroscopic redshift distributions. The distributions are 
in excellent agreement to z ~ 1.5. At 1.5 < z < 3, the pho- 



tometric redshifts are contaminated by low redshift galax- 
ies (see Section 5.2) and the paucity of spectral features 
in UV makes spectroscopic redshift measurement difficult 
(' Le Fevre et al. E 005al . Both effects explain the difference 
between the photometric and spectroscopic redshift distri- 
bution at z > 1.5. 

We are able to identify peaks in the photometric red- 
shift distribution which are clearly associated with peaks 
in spectroscopic redshift distribution. We smooth the 
spectroscopic and photometric redshift distributions us- 
ing a sliding window with a step Az — 0.2. The ratio 
between the observed redshift distribution obtained with 
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Fig. 11. Az as a function of redshift in four apparent mag- 
nitude bins (shown for clarity). We report the 3cr error 
bars on the photometric redshift estimate. 

a step Az = 0.01 and the smoothed redshift distribution 
shows three peaks at zp ^ 0.31,0.61,0.88 in the photo- 
metric redshift distribution, in excellent agreement with 
peaks identified at zs ^ 0.33, 0.60, 0.89 in the spectro- 
scopic redshift distribution. The significance of the peaks 
is lower by a factor two in the photometric redshift sample 
since the peaks are broadened by the uncertainties on the 
photometric redshift estimates. 

6. Added value of each multi-colour data set 

We now compare the reliability of photometric redshifts 
computed from either a BVRI dataset using the VVDS- 
CFH12K photometry or u*g'r'i'z' using the CFHTLS- 
MEGACAM, and quantify the useful range of photometric 
redshifts for each of these datasets. In addition, we explore 
the added value of different bands to the accuracy and re- 
liability of photometric redshifts. 

6.1. Added value of u* and z' bands 

Removing successively the u* and the z' bands, we inves- 
tigate possible systematic trends if these bands are not 
available or if they are shallower. 

Figure El (top left-hand panel) shows the photometric 
redshifts computed without u* band data. Only ~ 70% of 
the photometric redshifts at zs < 0.4 are recovered which 
should be compared with ^ 90% using the u* band. Since 
the filter system is no longer sensitive to the Lyman break, 
a large fraction of low redshift galaxies contaminates the 
zp > 3 redshift range. This test shows the importance of 



a deep u* band to constrain the photometric redshifts at 
z < 0.4 and z > 3. 

Figure lT^ (top right panel) shows the photometric red- 
shifts computed without z* band data. Most of the pho- 
tometric redshifts at zs > 1 are estimated at zp < 1. We 
observe an accumulation of photometric redshifts around 
zp ~ 0.8 — 0.9. This trend is expected since the filter sys- 
tem is no longer sensitive to the Balmer break at z > 1. 
In this case the use of the photometric redshifts is prob- 
lematic even at z < 1 without z' band data. 

6.2. Photometric redshifts from the VVDS imaging 
survey alone 

The VVDS multi-colour survey was carried out in B, V , 
R and / bands over 10 deg^ IjLe Fevre et al. E 004a') us- 
ing the CFII12K wide-field mosaic camera at CFHT. The 
accuracy of the photometric redshifts using only BVRI 
is presented in Figure 1131 (bottom right panel) . Since the 
VVDS photometric survey is shallower, the quality of the 
photometric redshifts is obviously worse than the results 
presented previously. 77 rises to 19.1% which is a factor 
of five greater than our best value. As we demonstrate in 
Section 6.3, the absence of deep u* and z' band data makes 
it difficult to compute photometric redshifts at zs < 0.4 
and zs > 1. However, even using only four broad bands, 
we recover 80% of the spectroscopic redshifts at Iab < 24 
with aiAz\/{i+zs) = 0.057. 

6.3. Photometric redshifts from the CFHTLS imaging 
data alone 

The deep CFHTLS survey consists in four fields im- 
aged over 3.2 deg^ in the u*g'r'i'z' filters. The qual- 
ity of the photometric redshifts computed using only the 
u*g'r'i'z' bands is displayed in Figure ITU ftop panel) for 
the CFHTLS-Dl. We find rj = 4.2% and cr|Az|/(i+zs) = 
0.040. These CFHTLS photometric redshifts are close to 
be as accurate as the photometric redshifts computed us- 
ing the full photometric dataset. Photometric redshifts 
for the other CFHTLS deep fields will be introduced in 
Section 7. 

6.4. The near-infrared sample 

Deep near-infrared observations in the J and K bands 
are available for a 160 arcmin^ l|Iovino et al. 2005) sub- 
area of the Dl field. This complete sub-sample of 3688 
galaxies at Kab < 23 represents a unique dataset in 
term of depth and area (it is one magnitude deeper and 
covers a three times larger area than the K20 survey, 
ICimatti et al. E 002V Near- infrared bands are crucial to 
constrain the photometric redshifts in the 'redshift desert' 
since the J band is sensitive to the Balmer break up to 
z 2.5 and enters in the K band at z > 3.8. The pho- 
tometric redshifts for galaxies selected at Kab < 23 are 
shown in Figure El (bottom left panel). We obtain the 
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Fig. 12. Comparison between the photometric redshift distributions and the VVDS spectroscopic redshift distributions 
on the CFHTLS-Dl field, for samples selected at < 23 (top left), < 24 (top right), Kab < 22 (bottom left) 
and Kab < 23 (bottom right). The black solid lines and the red dashed lines correspond respectively to the estimate 
with and without using the prior on the redshift distribution. These distributions are compared with the spectroscopic 
redshift distributions (shaded histograms) from the VVDS sample, originally selected at Iab < 24. To maintain the 
same vertical axis, the redshifts distributions are divided by a factor of two for i' < 24, K < 23. 



Magnitude cut 


I* < 23 


i* < 24 


Ks < 22 


Ks < 23 




0.76 


0.90 


0.90 


1.07 


% at 2 > 1 


13% 


28% 


28% 


43% 



Table 3. Median redshifts and fraction of galaxies at z > 1 for samples selected according to < 23, 24 and 
Kab < 22, 23 in the CFHTLS-Dl 
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Dl D2 D3 D4 



17.5 < i'AB < 22 0.49 0.48 0.49 0.48 

22 < i'^B < 23 0.75 0.75 0.75 0.81 

23 < i'AB < 24 0.89 0.86 0.86 0.95 

24 < i'^B < 25 1.00 0.95 0.94 1.03 

Table 4. Median redshifts in the four CFHTLS deep fields 
(columns) for samples selected according to 17.5 < < 
22, 22 < i'^B < 23, 23 < i'^^ < 24, 24 < i'^g < 25 from 
the top to the bottom, respectively. 



most reliable photometric redshifts on this sub-sample 
with 77 = 2.1% and <J\Az\/{i+zs) — 0.035. 

We present in Figure IT^ the redshift distributions for 
samples selected with Kab < 22 and 23. Comparing the 
photometric and spectroscopic redshift distributions, we 
see a large difference in the high redshift tail which is 
explained by the colour incompleteness caused by the 
Iab ^ 24 selection function of the spectroscopic sam- 
ple. The median redshifts and the fraction of galaxies 
with z > 1 are given in Table 13 As expected, we find 
that near-infrared selected samples are more efficient to 
target a high redshift population than i' selected sam- 
ple. We find 43% of the galaxies at z > 1 for a sample 
selected at Kab < 23. As previous K selected surveys 
IjCimatti et al. E 002. ISomerville et al~b 004') have found, a 
large population of galaxies at z > 1 is observed. 

7. Photometric redshifts in the CFHTLS "Deep 
Fields" Dl, D2, D3 and D4 

We finally use the photometric redshift calibration derived 
from the CFHTLS-Dl field and the VVDS spectroscopic 
sample to derive photometric redshifts for all fields of the 
CFHTLS deep survey. 

The three CFHTLS deep fields D2, D3, D4 have been 
imaged with the same instrument and are reduced ho- 
mogeneously in exactly the same way as the Dl field 
IjMcCracken et ar~b 006. in preparation). We therefore as- 
sume that we can measure the photometric redshifts for 
these fields in the same manner as we have for the Dl field. 
As a consistency check, we use 364 spectroscopic redshifts 
from the DEEPl survey ( [Phillips et al . 1997) which are 
in the D3 field. This allows us to test blindly the quality 
of these photometric redshifts without any additional cal- 
ibration. The comparison is shown in the bottom panel of 
Figure in We find 77 = 3.8% and cr|A^|/(i+^s) = 0.035 at 
^'ab — 24 and z < 1.5, without any systematic trend. We 
therefore conclude that our calibration method derived 
from Dl can be applied to the other CFHTLS deep fields. 

The photometric redshift quality using only the u* , g', 
r' , i' , z' has already been discussed in Section 6.3, but 
only for the CFHTLS-Dl field and for i'^g < 24. Since 
we have already demonstrated in Section 5.3 that the la 
error bars are representative of a measurement at 68% of 
confidence level, we use the la error bars to quantify the 



accuracy of the photometric redshifts in the different fields 
and beyond the spectroscopic Hmit. Figure [Tsl shows the 
fraction of photometric redshifts with a la error bar left 
than 0.15 x (1 -I- z). The best constraint is obtained on 
the CFHTLS-Dl field and gradually declines for the D4, 
D3 and D2. The constraint on the photometric redshifts 
is the lowest on the D2, which is expected since the total 
exposure times in the u* and z' bands are respectively 
7.7 and 1.7 times lower for the D2 field than for the the 
Dl field. We note that the specific trends described in 
Section 6.1 and shown in Figure [T^ could partially affect 
the photometric redshift estimates for the CFHTLS-D2 
given than this field has substantially shallower u* and z' 
data. The other significant trends observed in Figure ITSl 
are expected from our previous comparisons: 

— the redshift range 0.2 < z < 1.5 is the most suitable 
for the 4 fields which is expected since this redshift 
range is constrained by the set of filters used. 

— the accuracy of the photometric redshifts decreases to- 
ward fainter apparent magnitudes, faster at i'^^ > 24. 
For 0.2 < z < 1.5, the fraction of galaxies with 
(Tzp(68%) > 0.15 X (l + z) remains greater than ^ 80% 
at i'^g = 25 in the CFHTLS-Dl field. 

We show in Figure[T?)lthe redshift distributions for the 
four CFHTLS deep fields. As expected, the median red- 
shift increases for fainter samples (Table 01 rising from 
Zm ~ 0.45 at i'^g < 22 to z ~ 1 at 24 < i'^g < 25. The 
median redshifts are in good agreement between the four 
fields. The redshift distribution in the D4 is shifted at 
higher redshift. We observe significant variations of the 
redshift distribution between the four fields. Figure ITtI 
shows the ratio between the redshift distribution in each 
field and the redshift distribution averaged over the four 
fields, using a redshift step of Az = 0.1. This ratio shows 
that the difference in redshift distribution can reach at 
most a factor 1.6 in a redshift bin A = 0.1 (at z=0.25 be- 
tween the D2 and the D4 fields). The average dispersion in 
the interval 0.2 < z < 1.5 and in a redshift slice A = 0.1 
is ~ 15%. We conclude that the cosmic variance can be 
important for fields of ^--^0.8 deg^. We note however that 
the four fields do not reach exactly same depth in all fil- 
ters which could be responsible for some of the differences 
between the redshifts distributions. 

8. Conclusions 

Using the unique combination of the deep u* g'r'i' z' multi- 
band imaging data from the CFHTLS survey supple- 
mented by shallower BVRI data from the VVDS imaging 
survey (and also by J and K data on a smaller sub-area) 
and VVDS first epoch spectroscopic redshifts, we have 
been able to obtain very accurate photometric redshifts 
on the CFHTLS-Dl field. We reach crAz/(i+2) = 0.037 
at i' < 24 and 77 = 3.7% of catastrophic errors (defined 
strictly as Az > 0.15(1 + z)). For the bright sample se- 
lected at lAB < 22.5, we reach crAz/(i+z) — 0.030 and 
7? = 1.7%. 
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Fig. 13. Comparison between spectroscopic and pliotometric redsliifts for different combinations of filters. The pho- 
tometric redshifts in the top left and right panel are computed without using the deep u* jz' band respectively. The 

bottom left panel shows the photometric redshifts for a near-infrarod selected sample computed using _B, V ^ R, /, u* , 
g', r', i, z', J and K bands. The bottom right panel shows the photometric redshifts obtained using the B, V, R, I 
bands from the WDS survey. 



This accuracy has been achieved by calibrating our 
photometric redshifts using a large and deep spectro- 
scopic sample of 2867 galaxies. We have established a reli- 
able calibration method combining an iterative correction 
of photometric zero-points, template optimisation, and a 
Bayesian approach. This method removes some obvious 
systematic trends in the estimate and reduces by a factor 
2.3 the fraction of catastrophic errors. 

We have investigated in detail the quality of photo- 
metric redshifts as a function of spectral type, apparent 



magnitude and redshift based on the comparison with 
the VVDS spectroscopic redshifts. This step is crucial for 
forthcoming scientific analysis. As expected we find that 
the apparent magnitude is the key parameter: the frac- 
tion of catastrophic errors increases by a factor 12 and 
the rms by a factor 1.7 between 17.5 < < 21.5 and 
23.5 < i'j^g < 24. The reliability of the photometric red- 
shifts also depends on the spectral type: half of the catas- 
trophic errors arc galaxies which are best fitted by a star- 
burst template type. The evolution of ?7 as a function of 
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Fig. 14. Az as a function of redshift. The photometric redshifts are computed using the CFHTLS fiher set u* , g' , r' , 
i', z'. The top and bottom panels present the photometric redshifts obtained on the CFHTLS-Dl and CFHTLS-D3 
fields respectively. 



redshift shows that the redshift range the most reliable 
for forthcoming scientific analysis is 0.2 < zp < 1.5 for 
the complete population of galaxies. This range can be ex- 
tended into the "redshift desert" when near-infrared data 
are available (although currently only 6% of the field is 
covered) . 

We present i' band selected redshift distributions at 
i'AB — 23 and < 24 which are fully consistent with 
the redshift distributions derived from the VVDS spec- 
troscopic redshifts. We show the ability of our method 
to correctly recover the redshift distributions, including 
even the identification of the strongest density peaks. We 
show that a near-infrared selected sample is very eflacient 
for the selection of high redshift galaxies, with 40% of 
the sample at z > 1 for K < 23. This robust K selected 
sample will be used to investigate the evolution of the 
stellar mass function IjPozzetti et al. IB006. in prepara- 
tion) which is a crucial test of the hierarchical model 
(e.g. lEauffmann fc Chariot I L998. ICimatti et 30 002. 
ISomerville et al."b 004) of galaxy formation. 

Finally, we have applied our robust photometric 
redshifts measurement code on all four CFHTLS dee p 
fields I http://www. cfht. hawaii. edu/Science/CFHTLSl . 
We measure photometric redshifts for an uniquely large 
and deep sample of 522286 objects at i'j^^ < 25 on 
3.2 deg^. We assess the accuracy of these photometric 
redshifts beyond the spectroscopic limits and we present 
the redshift distributions in these four deep fields showing 



that cosmic variance effects are present at the 15% level 
for fields of size 0.8 deg^. 

All photometric redshifts and input photometric cata- 
logues are made publicly available. 
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Fig. 17. Ratio between the redshift distributions in each 
field and the redshift distribution averaged over the 4 
fields. Symbols for each field are the same than Figure [TBI 



